Origin of charge density wave in the coupled spin ladders in Sri^a-Ca^Ci^O/ii 
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We formulate a d-p multiband charge transfer model for CU2O5 coupled spin ladders, relevant 
for CU2O3 plane of Sri4_ a; Ca I Cu2404i, and solve it using Hartree-Fock approximation. The results 
explain that (i) the charge density wave (CDW) with its periodicity dependent on doping is stabilized 
by purely electronic many-body interactions in a single spin ladder and (ii) the inclusion of the 
interladder interactions favors (disfavors) the stability of the CDW with odd (even) periodicity, 
respectively. This stays in agreement with recent experimental results and suggests the structure of 
the minimal microscopic model which should be considered in future more sophisticated studies. 

PACS numbers: 74.72.-h, 71.10.Fd, 71.45.Lr, 75.10.Lp 
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Sr 1 4_ ;E Ca a ;Cu2404i (SCCO) is a layered material with 
two distinctly different types of two-dimensional (2D) 
copper oxide planes separated by Sr/Ca atoms>i (i) the 
planes with almost decoupled Cu02 chains and (ii) the 
CU2O3 planes formed by C112O5 coupled ladders (see Fig. 
[TJ). The latter ones exhibit the non-BCS superconducting 
(SC) phase for x = 13.6 in SCCO under pressure larger 
than 3 GPa^ or a spin-gapped insulating charge density 
wave (CDW) state in broad range of x and under normal 
pressured By means of the resonant soft x-ray scattering 
it was founds that this CDW is driven by many-body in- 
teractions (presumably just Coulomb on-site interactions 
since the long-range interactions are screened in copper 
oxides 5 ), and it cannot be explained by a conventional 
Peierls mechanism. Hence, the observed competition be- 
tween the CDW (also referred to as the "hole crystal" 
due to its electronic origin) and SC states in spin ladders 
resembles the one between stripes and the SC phase in 
CuC>2 planes of a high T c superconductor (HTS), which 
makes the problem of the origin of the CDW phase in 
SCCO both generic and of general interest. 

Furthermore, recently it has been found 5 that the only 
stable CDW states are with period A = 5 for x = 0, and 
with period A = 3 for x = 11 (and with a much smaller 
intensity for x = 10 and 12), while such a CDW order has 
not been observed for 1 < x < 5. These striking results, 
which contradict the previous suggestion 3 that the CDW 
order occurs in the entire range of < x < 10, need 
to be explained by considering hole density per Cu site 
increasing with x. Recently 7 a much higher hole density 
in the ladders was reported than believed before^ i.e., 
rih = 1-20 (number of holes/Cu ion) for x = 0, nh = 1.24 
for x = 4, and nu = 1.31 for x = 11. The aim of this 
paper is to explain theoretically these puzzling features 
of the CDW order using the above hole density. 

On the one hand, it is widely believed^ that a two-leg 
spin ladder described by the t-J model captures the es- 
sential physical properties of SCCO. The idea that merely 
on-site Coulomb interactions could lead to charge or- 
dering was already suggested by White et al^ using 
density matrix renormalization group (DMRG) — they 



found that a CDW of period A = 4 is the (possibly spin 
gapped) ground state at rih = 1.25. It is, however, quite 
remarkable that such a CDW has not been observed** On 
the other hand, the validity of the t-J model for CU2O5 
coupled spin ladders is not obvious since: (i) unlike the 
Cu02 plane of a HTS, a single CU2O5 ladder lacks the 
Dnh symmetry making the Zhang-Rice (ZR) derivation-^ 
of the t-J model questionable and (ii) CU2O5 spin lad- 
ders are coupled through the on-site Coulomb interac- 
tions between holes in different 0(2p) orbitals, so new 
interactions could arise. 

This suggests that the multiband charge transfer 
modeli 2 - adapted to the CU2O5 ladder geometry, similar 
to those introduced earlier for Cu02 planes^ 3 - or CUO3 
chains^ of HTSs, could be more appropriate to cap- 




FIG. 1: (Color online) Schematic picture of two coupled 
CU2O5 ladders (white and grey) with a CDW order of pe- 
riod: (a) A = 3 and (b) A = 4. The Cu 2 5 unit cell with 
two ?>d x i_ y 2, three 2p x , and two 2p y orbitals is indicated by 
dashed line. The arrows stand for hole spins in Cu and O 
orbitals, with their (large) small size corresponding to +1.0 
(+0.25) hole charge. The ovals show rungs with enhanced 
hole density in the CDW phase. The dotted ovals in the grey 
ladder of (b) show the two possible degenerate states, see text. 
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ture the essential physical phenomena. As parameters 
the charge transfer model includes: the energy for oxy- 
gen 2p orbital A, the d-p hopping t between the nearest 
neighbor Cu and O sites, and the on-site Coulomb repul- 
sion U (Up) on the Cu (O) sites, respectively. By solving 
this model in the Hartree-Fock (HF) approximation, we 
investigate whether: (i) the Coulomb on-site repulsion 
stabilizes the observed CDW of the respective period A 
for a given number of holes n/j, (ii) the different stability 
of the CDW order with odd/even period follows, and (iii) 
the ZR singlets are formed in the spin ladder geometry. 
We consider the charge transfer model in hole notation 
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where \t m j\ — \t m l\ = t, a = x,y, and a = —a for 
a =t, I- The parameter rj = Jh/U p = 0.2 stands for 
a realistic value of Hund's exchange (U p is the intraor- 
bital repulsion), 5 and e — 0.92 yields the correct orbital 
energy at bridge positions.— The model of Eq. {1} in- 
cludes seven orbitals per CU2O5 ladder unit cell (see Fig. 
[1]): two Cu(3d x 2_ y 2 = d) orbitals on the right or left 
(R or L) leg, two 0(2p y = y) orbitals on the R/L leg, 
two 0(2p x = x) side orbitals on the R/L leg, and one 
0(2p x = b) bridge orbital on the rung of the ladder. We 
emphasize that the last two terms in Eq. ([1]) account 
for interladder interaction - the holes within two differ- 
ent orbitals on a given oxygen ion in a leg belong to 
two neighboring ladders (shown as white/grey orbitals 
in Fig. [1]), and are described by charge operators nj X ( y ) a 
with/without tilde sign in Eq. |T]). This makes the model 
Eq. ([T]) implicitly 2D, though the band structure is one- 
dimensional (ID) when the interoxygen hopping t pp ' is 
neglected (in fact,— t pp i <C t). A priori, hj a should be 
treated as particle number operators belonging to the 
Hilbert subspace of the neighboring ladder, resulting in 
a 2D many-body problem. Here we simplify it and treat 
Pj a = (fija) as "external" classical fields adjusted to the 
symmetry of the CDW state, which are self-consistently 
determined within the HF approximation. 

We have solved the Hamiltonian ([1]) for various val- 
ues of the model parameters {U, A,U P }, and for three 
different hole densities 7 rih — 1.20, 1.25, 1.33 (which cor- 
respond to the actual filling in SCCO in the range of 
< x < 11) using HF approximation, i.e., we decou- 
ple n J>T n i/4 -> (n Jt2l )n jflL + n.j^(n ml ) - K >T )(n J> ;), 
where /1 — d,x, y, b. The ground state was found by 
diagonalizing the resulting one-particle Hamiltonian in 



real space for a single ladder with 60 unit cells, sepa- 
rately for spin up and spin down. The classical fields 
{Pja} and {(rij^a)} were determined self-consistently 
with the initial values for these fields suggested by recent 
experiment j 7 - see Fig. [T] While a uniform spin density 
wave (SDW) is stable for n^ — 1.0, one finds a CDW su- 
perimposed on the SDW order for realistic hole densities 
nil > 1.20. The stability of this composite order follows 
from the ID polaronic defects in the SDW state. We 
limit the present analysis to the stability of this particu- 
lar CDW phase, while we do not study here the possible 
competition with other phases^ 

For each state we evaluate the CDW order parameter 
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where A is the period of the CDW state, and the ZR 
" dispersion" defined with respect to the hole density dis- 
tribution for an "idealized" ZR singlet state (no = 0.25), 
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Here and in what follows by "rung" we mean the "rung 
with enhanced hole density" which consists of seven O 
(four y, two x and one b) orbitals and two Cu orbitals 
(see the ovals in Fig. [I]). Hence, in both above defini- 
tions the mean values of the particle number operators 
are calculated for these rungs (i £ rung) or for all re- 
maining sites (i rung). Note that in the ideal CDW 
phase (shown in Fig. [T]) p = 2 and a 2 = 0, irrespectively 
of the actual period A. We also introduce rung hole den- 
sities on O and Cu sites 

n p = Y, ( n ib + n ix + n iy ) , n d = ^ (nid)- (4) 

i£rung iGrung 

Similarly, magnetic order parameters are 
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where the magnetization for orbital /1 at site i is rn^ = 
( n in]- ~ n i^i)- We recall that when holes on the rungs 
form two localized ZR singlets next to each other, then 
n d = nid — 2, n p ~ 2, and m p ~ 1.5, see Fig. [1] 

First, we investigate the onset of the CDW phase in a 
single ladder of Fig. [T]by assuming U p — 0. In the charge 
transfer regime (for A = 3t following Ref. Il7t ) the CDW 
is stable already for U > t with periods: A = 5 for n^ — 
1.20, A = 4 for n h = 1.25, and A = 3 for n h = 1.33 [Fig. 
Ufa)] . For higher values of the on-site Coulomb repulsion 
U, p first increases quite fast irrespectively of the actual 
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FIG. 2: (Color online) Characterization of the CDW ground 
states obtained with U p = for increasing [/ (left, A = 3i) 
and A (right, U = 8£): (a), (b) CDW order parameter p, and 
(c), (d) ZR singlet dispersion a 2 , for A = 5,4,3 shown by 
solid, dashed, and dotted lines, respectively; (e), (f) for A = 5 
charge (magnetization) in the rung on Cu sites shown by solid 
(dotted) line and on O sites shown by dashed (dashed-dotted) 
line, see Eqs. I©-©. The realistic values (Ref. [T3) of U = 8t 
and A = 3t are marked by vertical lines. 

CDW period, and next saturates at p ~ 1, being only 
about 50% of the maximal value p = 2 (a weak decrease 
of p for U > 6t follows from the charge redistribution). 
In particular, such a CDW order is robust for the widely 
accepted value of U = 8t for copper oxide ladders^ 

In the strong coupling regime of U > At the CDW is 
formed by holes distributed as in the ZR singlets since 
then a 2 ~ 0.05 is indeed very small for all periods [Fig. 
[He)]. This is also visible in Fig. [5Je) where, in this 
regime, both the number of holes on O sites (n p ) and 
on Cu sites (n^) in the rungs are not far from their val- 
ues in the localized ZR states. Note that the minimum 
of u would correspond to n p — which further moti- 
vates the definition of Eq. (J3J). We can also probe the 
ZR character of holes forming the CDW by looking at 
the magnetization of holes in the rungs, cf. Fig. [DJe). 
The magnetization wid grows with increasing U and for 
large U ~ 12t it is still around 30% smaller than for lo- 
calized ZR singlets. However, even in this range of U 
the magnetization on the O sites m p is quite small and 
much below the value for ideal ZR singlets (around 70% 
smaller). This confirms that the subtle nature of the ZR 
singlets can be only partly captured within the HF ap- 
proach. 

Remarkably, changing the value of A for fixed U = 8t 
does not destabilize the CDW [Fig. IHb)] irrespectively 
of the period. This suggests that the charge order is 



FIG. 3: (Color online) Charge gap Aqdw in the CDW ground 
state obtained for increasing U, with period A = 5, 4, 3 (solid, 
dashed, and dotted line) . The inset shows the density of states 
N(E) for U = 8t and A = 5, with small CDW gap near the 
Fermi energy (dashed line). Parameters: A = 3t, U p = 0. 



triggered by the on-site Coulomb repulsion. However, 
the character of the holes forming the CDW changes and 
a 2 is small (a 2 ~ 0.07) only as long as A is large [Fig. 
[2^d)] . This is also visible in Fig. [^f) where a similar 
discussion as the one concerning Fig. [^e) applies. 

To gain a deeper understanding of the results we calcu- 
lated the charge gap Aqdw as a function of the Hubbard 
U, see Fig. [3] One finds that the gap decreases with the 
CDW period which explains the behavior of p observed 
for different periods [Fig. UK a)]. In general the depen- 
dence of Acdw on U qualitatively mimics the relation 
between p and U which suggests that the CDW gains 
stability when an insulating state is formed. Indeed, the 
electronic density of states N(E) (inset of Fig. [3} shows 
well developed lower and upper Hubbard bands (LHB 
and UHB) separated by an oxygen band, with a small 
CDW gap in the latter band. Altogether, one finds that: 
(i) the Coulomb interaction U can stabilize the CDW in 
the CU2O5 ladders, (ii) the CDW phase can be viewed 
as an equidistant distribution of the ZR singlet states in 
the relevant parameter regime, and (iii) all of the stable 
periods (even and odd) behave similarly. 

Next, we investigate the influence of the interlad- 
der coupling. At finite U p the "external" fields pj a = 
(nj+x,a) m Eq- © contribute and were self-consistently 
determined by iterating the HF equations. Thereby, the 
symmetry of the CDW state was chosen in such a way 
that the rungs were translated by A Cu-0 lattice con- 
stants (A odd) in the neighboring ladders to maximize the 
distance between them (Fig. [T]) , which minimizes the HF 
energy. For even A = 4 the numerical calculations per- 
formed with the realistic parameters^ for CU2O5 ladder 
(U = 8t and A = 3t) confirmed that two states shown by 
dotted ovals in Fig. [T]are degenerate, as expected. The 
effect of the interladder interaction U p was identified by 
comparing the ground states derived separately in two 
cases: (A) with pj a = 0, i.e., using only the (intraor- 
bital) repulsion between oxygen holes on the considered 
ladder; (B) by implementing the "external" fields {pj a } 
calculated self-consistently, i.e., including both the in- 
traorbital and interorbital Coulomb repulsion between 
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FIG. 4: (Color online) The CDW ground state for increasing 
U p : (a) CDW order parameter p and (b) ZR singlet dispersion 
a 2 , for A = 5, 4, 3 shown by solid, dashed, and dotted lines 
(squares, triangles, and circles) in case A (B), see text; (c) 
for A = 3 charge on Cu (O) sites in the rung shown by solid 
(dashed) line in case A and by squares (circles) in case B, see 
Eq. @; (d) for A = 4 charge in different y orbitals (riy, n y i, 
and n y 2, shown by diamonds, triangles down, and up) in the 
rung in casefB), see text. Vertical lines mark the realistic 
value (Ref. 0]]) of U p = St. Parameters: A = St, U = St. 



holes on oxygen sites. 

One finds that in case A the CDW order parameter p 
decreases in a similar way for all periods, cf. Fig. 0Ja), 
as well as for even period (A = 4) when the interladder 
coupling is switched on (case B). Remarkably, a quali- 
tatively distinct behavior is found for odd periods - here 
the interladder coupling supports the onset of the CDW 
phase and the order parameter either saturates or even 
increases with increasing strength of the on-site repulsion 
Up (as for A = 3), see also Fig. H£c). In fact, the inter- 
ladder coupling enhances the hole density in the rungs. 

Another striking effect is the qualitatively distinct be- 
havior of the ZR dispersion a 2 for odd and even periods, 
see Fig. IDJb). While for period A — 4 switching on the 
interladder coupling [B) drastically increases tr 2 with re- 
spect to the single ladder case (A), the results are pre- 
cisely opposite for odd periods A = 3, 5. Furthermore, 
this increase of a 2 with U p in case (B) is large for even 
period - its value ~ 0.1 found for large (but still realistic) 
Up ~ 3.5i is comparable to the value of the ZR dispersion 
for a single ladder with A ~ t [Fig. [21(d)] , where we do 
not expect stable ZR singlets. This large increase of a 2 



in this case follows from the geometrical frustration of 
the CDW state, as for even periods the two y orbitals in 
the same rung are not equivalent [one of them (say yl) 
is closer than the other one (say y2) to the rung in the 
neighboring ladder], as shown in Fig. BJd). We have also 
verified that the mean hole density rig — \{n y i + n y 2) 
almost does not change when the interladder coupling is 
switched off (not shown). 

Thus, we conclude that the interladder interaction: (i) 
supports the CDW states with odd periods A = 3, 5 and 
slightly disfavors the frustrated CDW state with even 
period A = 4, (ii) destabilizes the homogeneous ZR-like 
distribution of holes in the rungs for period A = 4. In 
contrast, experimentally one finds that in SCCO with 
x — 4 (rih ~ 1.25) the holes are distributed isotropically 
over O sites in the rung^ but the CDW is unstable^ We 
suggest that, since in reality the ZR singlets are much 
more rigid than the present classical ZR states (as the 
energy gain due to quantum fluctuations and phase co- 
herence are not captured in these states) and in reality 
the system is less prone to order than in the HF approxi- 
mation, the interladder interactions in the model Eq. ([T]) 
would indeed destabilize the CDW with even period. 

In summary, we have shown that the CDW combined 
with the SDW can be stabilized in the spin ladders of 
SCCO merely due to on-site Coulomb repulsion on Cu 
sites. The presented results explain the experimentally 
observed CDW states with odd periods for x = and 
x = 11, and provide a theoretical explanation why the 
CDW states with even period could not be observed^ In 
addition our results suggest that an extension to the two 
leg ladder t-J model used in Ref. Qj] is needed to capture 
the subtle properties of the CDW states in SCCO. The 
simplest extension might be to consider a pair of ladders 
in a t-J model plus an interladder Coulomb repulsion 
to represent the physics described in our study. We are 
looking forward to future studies of this kind using the 
much more sophisticated DMRG-like methods. 

This work was supported by the Polish Ministry of Sci- 
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Note added in proof: In contrast to the present charge 
transfer model, Roux et al^ have found that the CDW 
state is unstable for other hole densities than rih = 1-25 
and rih = 1-50 in their recent DMRG and bosonization 
studies of the t-J model for the single ladder with 1.0 < 
n h < 1.5. 
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